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Qh| Abstract 

OhI We calculate the rate of neutron-mirror-neutron oscillations for 

ultracold neutrons trapped in a storage vessel. Recent experimental 
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bounds on the oscillation time are discussed. 



> 

C^ : 1 Introduction 

in 

in 

O ' During the last couple of years we are witnessing the revival of the interest to 

t:;J^ . the "mirror particles", "mirror matter" and "mirror world". The idea of the 

O I existence of the hypothetical hidden sector to compensate mirror asymmetry 

Q I was first explicitly formulated in [1]. The subject has a rich history - see the 

review paper [2|. The present wave of interest to mirror particles has been to 
a great extent initiated by the quest for neutron-mirror-neutron oscillations 
{n—n'). It was conjectured that n—n' oscillations may play an important role 
in the propagation of ultra high energy cosmic rays and that the oscillation 
time Tosc may be as small as Tosc ~ Is [3] . Implications of mirror particles for 
cosmology and astrophysics were discussed in a number of papers e.g., [1]. 
Last year the first experimental data on n — n' transitions were published 
with the results Tqsc > 103s [5j and Tosc > 414s [6]. Possible laboratory 
experiments to search for n — n' oscillations were discussed in [7j. 

Experimental results [5l [6] were obtained using the ultracold neutrons 
(UCN), i.e., neutrons with the energy E < 10""^ eV storaged in a trap. 
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Previously similar experimental setup was used in the search for neutron-anti- 
neutron oscillations (see [S] and references therein). The crucial difference 
between n — n' and n — n oscillations is that n' freely escapes from the 
trap while fi either annihilates on the trap walls or gets reflected. Therefore 
formalism developed for the n — n oscillations cannot be adjusted to treat 
n — n' transitions. Still the two processes have a common point. This is 
a problem of a correct quantum mechanical description of the UCN wave 
function (w.f.). Most often it is assumed that the w.f. of the bottled UCN 
corresponds to a stationary state of a particle inside a potential well [9l [10] . 
Alternatively, other authors |llj describe oscillations of the trapped neutrons 
in the basis of the free plane waves. Both pictures do not correspond to the 
physics of real experiments. The process proceeds in time in three stages. 

At the first stage the filling of the trap takes place, then the neutrons 
are kept inside the trap during the storage time (hundreds of seconds), and 
finally neutrons leave the trap to the detectors. Therefore the w.f. undergoes 
a complicated evolution which hardly can be described without resorting to 
approximations. We shall first evaluate the neutron-mirror-neutron oscilla- 
tions using a stationary w.f. as the initial state w.f. Then we shall do the 
same using wave packet instead of a stationary w.f. 

The paper is organized as follows. We start in Section 2 with the analysis 
of the oscillations in the stationary w.f. approach. Transitions take place 
from one of the trap eigenstates. In Section 3 transitions are considered 
in presence of a superimposed magnetic field. A general equation for the 
transition rate is derived and the limits of weak and strong field are con- 
sidered. Section 4 is devoted to the wave packet formalism. The evolution 
of the UCN wave packet (w.p.) is encoded using the trap Green's function. 
Neutron-mirror-neutron transition rate is calculated. In Section 5 the main 
conclusions are presented and open problems are formulated. Appendix con- 
tains comparison between the infinite and finite well models. 

2 Stationary Wave Function Approach 

The problem of neutron-mirror-neutron oscillations in free space can be 
solved by diagonalization of the time-dependent two-channel Schrodinger 
equation with the result [12] 



\^nW = ^^j^exp(-r^t) sin^ ^^VIJ^T^t^ , (1) 



where uj = En — En' = \fJ'n\B is the energy difference between neutron and 

mirror neutron due to superimposed magnetic field (mirror neutron does 

not feel "our" magnetic field), e = t~^1 is the mixing parameter, Tp is the 

neutron /3-decay width. In arriving at ([1]) the spatial part of the w.f. was 

factored out making use of the fact that in free space the w.f.-s of n and n' 

are of the same form. In the trap, however, the situation is different: the 

neutron is confined while for the mirror neutron the trap walls do not exist. 

As already mentioned in the Introduction, the description of the trapped 

UCN is a nontrivial problem. The naive guess would be that inside the trap 

the neutron w.f. corresponds to a discrete eigenstate. Here we assume that 

the neutron w.f. is that of a particle in a potential well with the boundary 

conditions corresponding (in the first approximation) to a complete reflection. 

In order to make calculations tracktable and transparent we shall consider 

the following simple model of a trap. Let it be a one-dimensional square well 

of width L = Im with walls at a; = and x = L, i.e., the potential of the 

form 

( V, x<0 

U{x) = <^ 0, < X < L . (2) 

[ V, x> L 

The height of the potential well depends on the material with the typical 
value V = 2 ■ 10^^ eV which will be used in what follows. For such a well 
the limit for stored UCN velocity is 6.2m/ s. The number of discrete levels 
in such a trap may be estimated as 



,^ LV2mV 10^ 

M ~ ~ . (3) 
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We choose the UCN energy to be ii^ = 0.8 ■ 10^^ eV. This energy corre- 
sponds to a level with quantum number j — 2- 10^. Positions and eigenfunc- 
tions of such highly excited states in a finite-depth potential are very close 
to the same quantities in the infinite well (except for the levels close to the 
upper edge of the well; we not consider such levels). The finite-depth cor- 
rections are considered in the Appendix. The eigenvalues and eigenf unctions 
for the infinite well are 

(Pj{x) = \ — sin kjX. (5) 



Another important quantity characterizing highly excited states is the 
classical frequency uici [H] 

TT^ , 27r ^^ _. 

where t^ is the time of the classical period and 

5Ej = Ej+i - Ej ^ 0.8 ■ 10"^^ eV 

is the level spacing. Levels with j ^ 1 are almost equidistant. In the 
semiclasical limit we may also define the trap crossing time r 

Tri niL , ^ 

r = ^ = -— ~ 0.26s (7) 

2 kj 

for j = 2 ■ 10^. Next we calculate the rate of {n — n') oscillations for the 
neutron at the j — th discrete level. The neutron and mirror neutron w.f.-s 
in a two-component basis are 

'fjix) = y-^sinfcjx i^j=ifj{x)i^j, (8) 

fvi^) = -^e^'"" ( ? ) = /p(^) (l), (9) 



/27r 
where — oo < p < +oo. The {n — n') system is described by the Hamiltonian 

The states ([8]) and ([9]) are the eigenstates of i^o, therefore it is convenient 
to use the interaction representation. The probability to find at time t a 
mirror neutron instead of a neutron reads 

Pnw = j dp\{~U\e^^^-i f^dt'Wint{t')Y(p,)\\ (11) 

where Wint{t) = e'^^"^We~'^^"^. In the first order of perturbation theory we 
get 

Pnw = Jdp\{fp\ J^ dt'W,nt{t'M)\^ = 



e'jdp\f^dt'e~^'^^^-^^^''W,\^,)\\ (12) 



fc2 2 

where Ej = ^, Ep = ■^. The time-dependent integral is a standard one 



rt 4 sin ^ " ^' 

u.(E,) = I / Me-'i-'--'^'r = ,^Lg.)J - (13) 



{Ep—Ej)t 



/o (Ep - E^f 

The overlap of the w.f.-s reads 

sAp) - \(fM)\' - ,^(pf! ,.). ^in^ (g^) ■ J = 1.2.... (14) 
From ([I2D, (USD and ([HD we obtain 

/ + 00 
dpgj{p)w{Ep). (15) 

It is convenient to change integration from dp to dEp taking into account 
that g{p) = g{-p). 
Then 

-2 /"jTp 9{Ep)w{Ep) 



Pnn' = 2me' / dEp ''' ''' ' ''^ (16) 

J p 

where the factor 2 comes from the fact that two plane waves e^*^^ correspond 
to the same energy Ep. Both functions g{Ep) and w{Ep) are peaked at 
Ep = Ej. According to ([HD and ([T3D the widths AE^ and AE^ of the 
corresponding maxima are 

AE^ ~ vr/r, AE^ ~ 47r/t, (17) 

with r being the trap crossing time. At times t ^ r we may substitute 
g{Ep)/p by its value at p = kj and take it out of the integral (iT6l) . From ^\M 
one gets g{Ej) = L/Att. The remaining integration in flT6|) can be extended 
to (— cxD < Ep < +cxd) yielding 27rt. Collecting all pieces together we obtain 

Pnn' = e^Tt. (18) 

At very short times t <^ t the function w{Ep) becomes smoother than 
g{Ep). Hence w{Ep) can be taken out of the integral fITBl) . The remaining 
integral is time- independent while w{Ep) r^ t'^. As a result Pnn' ~ ^^^^ and 
we can not define the transition probability per unit time [12] . On the other 
hand, Eg. (1181) is valid only for times shorter than the neutron j3 -decay time tfs 
since we have define the eigenstate ([8]) neglecting the /3-decay. The condition 
T <^t <^tp was with a fair accuracy satisfied in experiments [5], E] . 



3 Stationary approach with magnetic field in- 
cluded 

The search for n — n' oscillations is experiments with bottled UCN is based 
on the comparison of UCN storage with and without superimposed magnetic 
field [3 E] ■ It is assumed that there is no mirror magnetic field in the labo- 
ratory and therefore the interaction of the neutron with magnetic field lifts 
the degeneracy and thus suppresses the oscillations. 

In magnetic field B the energy of the trapped neutron becomes equal to 

£, = |- + ,B, (19) 

where /i = — /i„ = 1.^1^in{,1^n = e/2mp). 

Inclusion of the magnetic field does not alter the functions w{Ep) and 
g{p) given by Eq.-s flT3|) and flT^ . There is, however, an important difference 
between our present considerations and the previous section. As we see from 



(TT9|) w{Ep) now peaks at p = ^Jk'j + 2m^B while the maximum of g{p) is 
as before at p = ^kj. As a result instead of flT6|) we obtain 



cos^Mi /i+2^, J = 1,3,... 



Pr,n' = { , , \ '— (20) 

(^5)V,/rT^ sin^^./rT^, J =2,4... 
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This equation can be simplified taking into account that the quantities 
inB) and k'^/2m differ by many orders of magnihides. In experiments [51 E] 
the value of the magnetic field varied in the interval (1— 2)nT < B < (few)/xT 
which corresponds to 10~^^ eV < fiB < 10^^^ eV, while /c|/2m ^ 10~^ eV 
(note that the unshielded Earth magnetic field corresponds to fiB ~ 3 • 10^^^ 

eV« 1^). 

Therefore Eq. fl20l) easily reduces to 

r. t sin^ (i/i-Br) , , 

where r is the trap crossing time. For our model of the trap described in 
section 2 we have r/2 ~ 2-10^^ eV~^. Therefore in the limit of weak magnetic 
field S ~ nT Eq. ([HD yields 

Pnn' ^ e^Tt, (22) 



as expected (see (IT5I) ). In the opposite limit of strong magnetic field B ~ 
(few) nT we have to take into account that the quantities r and B in (1211) ex- 
perience fiuctuations leading to rapid oscillations of the function sin^ (|/ii?r). 
In particular, the crossing time r may vary either due to changes of L at 
each cossision, or due to variations of the neutron velocity. Substituting 
the rapidly oscillating quantity in fl?I]) by its mean value equal to 1/2 we ob- 
tain the equation describing the neutron-mirror-neutron transitions in strong 

magnetic field 

2t 

Pnn' = ^"j^^B)H- (23) 

4 The wave packet approach 

We now turn to the question formulated in the Introduction, namely to the 
problem of the UCN w.f . evolution and to the calculation of the oscillations in 
the wave packet approach. In order to get physically transparent results and 
to avoid numerical calculations suited to a concrete experiment we assume 
that UCN coming to the trap from the source are described by the Gaussian 
wave packets (w.p.) [13]. 

The w.p. moving from the left and for t = centered at a; = is given 
by the following expression 

v^,(x, t = 0) = (7ra2)-V4exp |-i^Z_^ + ^kx\ , (24) 

where a is the width of the w.p. and k is its central momentum. The 
normalization of the w.p. fl23l) corresponds to one particle in the entire one- 
dimensional space, 

r+oo 

/ t/x|^fc(x,t = 0)|^ = 1. (25) 

Let the UCN energy be equal to the value chosen in Section 2, E = 
0.8 • 10^''' eV, and let the beam resolution be equal to AE/E = 10"'^. Thus 
the set of parameters to be used i^j 

E = 0.8- 10"^ eV, A = -^ ~ 10"^ cm, a ~ 3.2 ■ 10"^ cm. (26) 

The condition a ^ X ensures the localization of the w.p. We remind that 
the above value of E corresponds to the level Ej with a very high quantum 



"'^The problem of the choice of the w.p. parameters will be addressed in the next Section. 



number j — 2 ■ 10^. Next we estimate the number of levels within AE". One 

AE viAk) L , 

Aj = = ^ ^ = — ~ 10^ 27) 

ujci ujci vra 

The large number of levels forming the w.p. is a necessary condition for 
the trapped w.p. to be localized (in free space this condition reads a ^ A, 
see above). The time evolution of the initial w.p. flM|) proceeds according to 
the following law 

^^(x, t)= j dx'G{x, t- x\ 0)^fc(x', 0), (28) 

where G{x,t;x',t') is the trap Green's function. In the infinite well approxi- 
mation we may use the spectral decomposition of the Green's function over 
the set of eigenfunctions ([5]) and write 

*fc(a:,t) = ^e-^^^V,(a;) / dx'ip*{x')^k{x' ,0). (29) 

The width of the w.p. (!29l) increases with time according to 



a' = a 



1 + ' * 



\ma? 



211/2 



( ^ 



\ 9 / ' 



where for our model the spreading time is ma"^ ^ 1.7 ■ 10~^s and t/ma^ ~ 
60t(s). A so-called collapse time tc [T^ corresponds to a' = L and is equal 
to tc — 500s. At t = tc the w.p. spreads uniformly over the entire well and 
the stationary regime considered in Section 2 sets in [15] . We note in passing 
that there is another time scale in the problem, the so-called revival time 
trev = 4mL^/7r ~ 2.10'''s when the w.p. regains its initial shape -see [H] and 
references therein. 

The initial w.p \E'fc(x, 0) contains only right running wave - see (l24l) . The 
trapped w.p. (129|) contains both right and left running waves, i.e., it correctly 
describes reflections from the trap walls. We assume that the point xq (see 
( IMl) is not in the immediate vicinity of the trap walls, i.e., xq is at least 
few times of a away from the walls. Then the integration in ( l29ll may be 
extended to the entire one-dimensional space . This yields 



/ + 00 / 

dx' ip*{x')'^ k{x' ,0) = i 




X < exp 



Qj \K fX/j 



+ i{k — kj)xo 



exp[...kj -^-kj...] 



(30) 



Then we can calculate the transition probability P„„/ following the pro- 
cedure described in Section 2. Instead of the w.f. (18]) we now have 






(31) 



with Fj{k) being the shorthand notation for the function Fj{k,kj; L,a,xo) 
defined by (130|) . The normalization condition for Fj{k) reads 



Ei^.(^)r = i- 



(32) 



In line with ( TTSl) and following the arguments presented after ( fT6l) we 
write 



Pnn'=e'T.^Ak)Fi*{k)e^^ 



HEi-Ej)t 



X 



X 



dp 



2 sin 



[Ep—Ej)t 



{Ep - Ej] 



2 sin 



{Ep-Ei)t 



{Ep - El 



Uv\v,){^i\fp)- (33) 



>(i) 



Consider first the contribution P„J, of the diagonal terms with j = I. We 
have 

with (r) being the weighted crossing time 



'3J^ 



(34) 



(35) 



and T(kj) = mL/kj. Next we turn to the contribution P^J„ of the nondiago- 
nal terms in (32). In this case we are dealing with a two-hump function with 
maxima at Ep = Ej and Ep = Ei. Therefore we may write 



p5 = 4.™= k^E^n*)^'"^'-"'" '"" 

Replacing summation over / by integration we obtain 



(Ej-Ei)t 



{E, - El 









e'{r') 



X 



(36) 



(37) 



where 



(r^) = Ei^.r-^ 



(38) 



Collecting the two contributions (1M1) and ( 1371) together we get the final 
result 



Pnn'=e\T)t[l + 






(39) 



5 Conclusions 

We have calculated the rate of neutron-mirror-neutron oscillations for trapped 
UCN. Two types of the UCN w.f.-s were used: the stationary solution for 
a particle inside a potential well and the Gaussian w.p. Calculations were 
performed in the first-order perturbation theory. This approximation is le- 
gitimate provided P„„' ^ 1. From (1221) and (1231) it follows that this condition 
holds for Tnn' ^ 7s and r„„/ ^ 0.1s in the weak and strong magnetic fields 
correspondingly. Obviously the first order perturbation theory describes the 
transition of the neutron into mirror neutron. The inverse process appears 
only in the second order in e. For the analysis of the experiments [Sj [6] first 
order perturbation theory is a fair approximation. 

The above analysis has been performed for a simple one-dimensional trap. 
We think that such a model correctly describes the principal features of the 
process. Generalization to the three-dimensional rectangular trap is triv- 
ial. The simplest way to generalize our result to the trap with arbitrary 
geometry is to substitute the crossing time r by the effective crossing time 
corresponding to a given trap geometry. 
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Experimental data [H E] were analyzed using the free space equation ([T]) 
with the time t being limited by the crossing time r. 

Equation ([1]) contains only time dependence since the spatial parts of n 
and n' w.f.-s were factored out using the fact that in free space the coordinate 
w.f.-s of n and n' have the same form. For bottled UCN the situation is dif- 
ferent. Neutron is confined inside the trap while mirror neutron freely crosses 
the trap walls. Therefore the use of the ([T]) to describe oscillations of trapped 
UCN seems questionable!. However our accurate approach justifies the anal- 
ysis of the experimental data based on Eq. ([I]) O [6] . This can be explained 
by the semiclassical character of the UCN motion inside the trap of macro- 
scopic size. In the stationary approach the typical UCN energy corresponds 
to the states with j 3> 1, i.e., to the semiclassical part of the spectrum. In 
the w.p. formalism the classical limit corresponds to (Aj) ~ {jY -^ cxo 
[16j . The w.p. flM|) . fl26l) is close to this limit. For UCN with extremely low 
energy, E < 10~^^ eV, the oscillation pattern changes. We shall consider this 
question elsewhere. The fraction of UCN with such energies in experiments 
is negligible. Another point which deserves a dedicated study is decoherence 
of the UCN w.f. and subsequent randomization of the oscillation process. 
This might occur due to dephasing of the w.f. caused by collisions with the 
trap walls and with the residual gas. 
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olation at Berkeley in September 2007. Also B.K. thanks V.A.Gordeev and 
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presentation of this work was made. Financial support from grants RFBR 06- 
02-17012, NS-4961.2008.2, NS-4568.2008.2, RFBR 07-02-00830 and RFBR- 
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Foundation for the financial support. 



^To describe the free space experiments by Eq.([T]) one still has to supplement it by a 
boundary condition at x = where the reactor is placed. Otherwise at i = 7rT„„'/2 the re- 
actor becomes a source of mirror neutrons. We are grateful to L.B.Okun and M.I.Vysotsky 
for drawing our attention to this. 
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6 Appendix 

Calculations presented above were performed for the infinite well model of 
a trap. Here we consider the finite potential and show that there is only a 
minor difference between the two models. Consider the potential well defined 
by Eq. ([2]). Matching the logarithmic derivatives of the w.f.-s at x = and 
X = L we obtain the eigenvalue equation 



k'jL = TTJ — 2 arcsin ^===, (A.l) 



(the notation kj is kept for kj = vrj/L). The small parameter in the problem 

is 

2 \i/2 



Expanding (]A.1|) with respect to 6 we obtain 



2 • 10"^ (A.2) 



Therefore the levels in the finite well are shifted relative to the infinite well 
levels by 

Ej -E'jC^A- 10"^^ eV. (A.4) 

From flA.3l) it follows that the spectrum in the finite well ([2]) is the same as 
in somewhat wider infinite well 

L' = L{l + 6). (A.5) 

In the finite well the w.f. penetrates into classically forbidden regions 
inside the trap walls. However neutron-mirror-neutron transitions inside the 
walls may be neglected since both the penetration depth d and the collision 
time Tcoii are small: d ~ 10~^ cm, Tcou ~ 10~^s. 
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